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Abstract 

We present a generalization of the resonant neutrino conversion in matter, including a 
random component in the matter density profile. The study is focused on the effect of such 
matter perturbations upon both large and small mixing angle MSW solutions to the solar 
neutrino problem. This is carried out both for the active-active v e — > v^ T as well as active- 
sterile v e — > v s conversion channels. We find that the small mixing MSW solution is much more 
stable (especially in 8m 2 ) than the large mixing solution. Future solar neutrino experiments, 
such as Borexino, could probe solar matter density noise at the few percent level. 



"Invited talk presented by A. Rossi at 17th Int. Conf. on Neutrino Physics and 
Astrophysics , Helsinki, Finland, 13-20 June 1996. To appear in the Proceedings. 



1. The comparison among the present experimental results on the observation of the solar 
neutrinos strongly points to a deficit of neutrino flux (dubbed the Solar Neutrino Problem (SNP)). 
The most recent averaged data jjj of the chlorine, gallium [] and Kamiokande experiments are: 

= (2.55 ± 0.25)SNU, R%» = (74 ± 8)SNU, R%* = (0.44 ± 0.06)i?f £ 95 (1) 

where R^a 95 1S prediction according to the most recent Standard Solar Model (SSM)by Bahcall- 
Pinsonneault (BP95)|] . 

It is now understood that the SNP cannot be explained through astrophysical/nuclear solutions 
||, [|]. From the particle physics point of view, however, the resonant neutrino conversion (the 
Mikheyev-Smirnov-Wolfenstein (MSW) effect) |J seems to explain successfully the present experi- 
mental situation || ^, ||, 

This talk deals with the stability of the MSW solution with respect to the possible presence of 
random perturbations in the solar matter density ||. 

We remind that in Ref.[^] the effect of periodic matter density perturbations added to a mean 
matter density p upon resonant neutrino conversion was investigated. There are also a number of 
papers which address similar effects by different approaches [11, 12]. 



Here we consider the effect of random matter density perturbations Sp(r), characterised by an 
arbitrary wave number k, 

5p(r)= [dk5p(k)srakr, (2) 



Moreover, as in Ref. |12f| , we assume that the perturbation 5p has Gaussian distribution with the 
spatial correlation function (£ 2 ) defined as 



(5p(r 1 )5p(r 2 )) = 2p 2 (e)L S(n-r 2 ), (f)^^l. (3) 



The correlation length Lq obeys the following relation: 

I free « I « Am (4) 

where Zf ree ~ 10 cm is the mean free path of the electrons in the solar medium and X m is the neutrino 
matter wave length. For the sake of discussion, in the following we choose to adjust Lq as follows: 

L = 0.1xA m . (5) 

The SSM in itself cannot account for the existence of density perturbations, since it is based on 
hydrostatic evolution equations. On the other hand, the present helioseismology observations cannot 
exclude the existence of few percent level of matter density fluctuations. Therefore, in what follow 
we assume, on phenomenological grounds, such levels for £, up to 8%. 

Before generalizing the MSW scenario, accounting for the presence in the interior of the sun of 
such matter density fluctuations, first we give a quick reminder to the main features of the MSW 
effect. 

2. The resonant conversion of neutrinos in a matter background is due to the coherent neutrino 
scattering off matter constituents |J. This determines an effective matter potential V for neutrinos. 
In the rest frame of the unpolarized matter, the potential is given, in the framework of the Standard 
Model, by 

V = (6) 

m p 

where Gf is the Fermi constant and Y is a number which depends on the neutrino type and on the 
chemical content of the medium. More precisely, Y = Y e — hY n for the v e state, Y = —\Y n for 

x For the gallium result we have taken the weighted average of GALLEX -R^f = (77 ± 
8 ± 5)SNU and SAGE = (69 ± 11 ± 6)SNU data. 



and u T and Y = for the sterile v s state, where Y e ^ n denotes the electron and neutron number 
per nucleon. For the matter density p, one usually consider the smooth distribution, as given by the 
SSM@,p,[l|. 



For given mass difference 5m 2 and neutrino mixing 8 in vacuum, the neutrinos z/ e 's, created in 
the inner region of the sun, where the p distribution is maximal, can be completely converted into 
u y (y = fi, t or s), while travelling to the solar surface. 
This requires two conditions ||: 

1) - the resonance condition. Neutrinos of given energy E experience the resonance if the energy 
splitting in the vacuum Sm 2 cos 28 /2E is compensated by the effective matter potential difference 
AV ey = V e — V y . It is helpful to define the following dynamical factor A ey 

A ey (r)= 1 -[AV ey (r)- 6 -^co S 29] (7) 

which vanishes at the resonance, A ey = 0. This condition determines the value 

Pres = (m p cos29/2\/2GF)(Y e — Y y )5m 2 j E which, in turn, implies a resonance layer Ar. 

2) - The adiabatic condition. At the resonance layer, the neutrino conversion v e — > v y is efficient 
if the propagation is adiabatic. This can be nicely expressed requiring the neutrino wavelenght A m 
to be smaller than Ar ]5[, 

, . , , 6m 2 sin 2 28 Rn , , 

a r = Ar/(A m ) res = — — r^>l, R ~ O.LR , 8 

4irE cos 28 

X m = — — , Ar = 2p res tan 28\dp/dr\~ 1 . 

A 2 ey + (5m 2 ) 2 sin 2 20/(16.6?) 



3. Now we re- formulate the neutrino evolution equation accounting for a fluctuation term 5p 



P 

component L\V ey ^ for the matter potential. The evolution for the v e — v y system is governed by 



superimposed to the main profile p. The perturbation level £ = -ff induces a corresponding random 



d_f v e \ = H e H ey Wi/ e ] 
l dt \ v y ) \ H ey H y )\v y )> [y) 

where the entries of the Hamiltonian matrix are given by 

5tn 2 ~ i 

H e = 2[A ey {t)+A ey {t)l H y = 0, H ey = —sm28, A ey (t) = -AV ey (t)£ . (10) 

Here the matter potentials read as: 

AV eKT) (t) = ^^ P (t)(l - Y n ) , AV es (t) = ^^p(t)(l - |y n ) (11) 

for the v e — > v^t an d v e — * v s conversions, respectively. (The neutral matter relation Y e = 1 — Y n 
has been used.) 

The system (||) has to be rewritten averaging over the random density distribution, taking into 
account that for the random component we have: 

<i 2 £ +1 > = 0, (A^A^h)) = K6(t - h), K (t) = (A 2 ey (t))L = ^AV 2 y (t)(e)L . (12) 

We have obtained (see || for more details) the following system: 

V{t) = 2H ey l{U) 

li(t) = -2A ey (t)l(U) - €k(\J)K(\J) 

±(t) = 2A ey (t)K(U)-£K(U)l(U)-H^(£V(U)-oo), (13) 



where V{U) = (|i^| 6 ), TZ(U) = (Re(u^)) and J(U) = (Im(z/fi/?)). Now the " dynamics " is 
governed by one more quantity i.e. the noise parameter k, besides the factor A ey . The quantity 
k can be given the meaning of energy quantum associated with the matter density perturbation. 
However, let us note that the MSW resonance condition, i.e. A ey (t) = remains unchanged, due to 
the random nature of the matter perturbations. The comparison between the noise parameter n in 
Eq. (|l^) and A ey (t) shows that K(t) < A ey (t), for £ < few %, except at the resonance region. As a 
result, the density perturbation can have its maximal effect just at the resonance. Furthermore, one 
can find the analogous of condition 2) (see Eq. (|8[) for the noise to give rise to sizeable effects. Since 
the noise term gives rise to a damping term in the system (|l3|), it follows that the corresponding 
noise length scale 1/k be much smaller than the thickness of the resonance layer Ar. In other words, 
the following adiabaticity condition 

a r = Ar K res > 1 , a r K, a r - — ^—^ . (14) 

tan z 20 

is also necessary. For the range of parameters we are considering, £ ~ 10~ 2 and tan 2 20 > — 10~ 2 , 
and due to the r.h.s of (Q), there results a r < a r . This relation can be rewritten as K res < 5H res , 
where 5H res is the level splitting between the energies of the neutrino mass eigenstates at resonance. 
This shows that the noise energy quantum is unable to "excite" the system, causing the level crossing 
(even at the resonance) flQ[| . In other words, it never violates the MSW adiabaticity condition. From 
Eq. (|l4| ) it follows also that, in the adiabatic regime a r > 1, the smaller the mixing angle value the 
larger the effect of the noise. Finally, as already noted above, the MSW non-adiabaticity a r < 1 is 
always transmitted to a r < 1. As a result, under our assumptions the fluctuations are expected to 
be ineffective in the non-adiabatic MSW regime. 

4. All this preliminary discussion is illustrated in the Fig. 1. For definiteness we take BP95 
SSM H as reference model. We plot V as a function of E/5m? for different values of the noise 
parameter £. For comparison, the standard MSW case £ = is also shown (lower solid curve). One 
can see that in both cases of small and large mixing (Fig. la and Fig. lb, respectively), the effect of 
the matter density noise is to raise the bottom of the pit (see dotted and dashed curves). In other 
words, the noise weakens the MSW suppression in the adiabatic-resonant regime, whereas its effect 
is negligible in the non-adiabatic region. The relative increase of the survival probability V is larger 
for the case of small mixing (Fig. la) as already guessed on the basis of Eq. (14). We have also 



drawn pictorially (solid vertical line) the position, in the V profile, where 7 Be neutrinos fall in for 
the relevant 5m 2 ~ 10 -5 eV 2 , to visualize that these intermediate energy neutrinos are the ones most 
likely to be affected by the matter noise. 

5. Let us analyse the possible impact of this scenario in the determination of solar neutrino 
parameters from the experimental data. For that we have performed the standard x 2 fit in the 
(sin 2 26, 5m 2 ) parameter space. The results of the fitting are shown in Fig. 2 where the 90% confidence 
level (C.L.) areas are drawn for different values of £ . Fig. 2a and Fig. 2b refer to the cases of v e —* z/„ )T 
and v e — ► v s conversion, respectively. One can observe that the small-mixing region is almost stable, 
with a slight shift down of 5m 2 values and a slight shift of sin 2 29 towards larger values. The large 
mixing area is also pretty stable, exhibiting the tendency to shift to smaller 5m 2 and sin 2 20. The 
smaller 5m 2 values compensate for the weakening of the MSW suppression due to the presence 
of matter noise, so that a larger portion of the neutrino energy spectrum can be converted. The 
presence of the matter density noise makes the data fit a little poorer: Xmin = 0.1 for £ = 0, it 
becomes Xmin = 0-8 f° r £ = 4% and even Xmin = 2 for £ =8% for the u e —> v^ T transition. 

The same holds in the case of transition into a sterile state (Fig. 2b): Xmin = 1 for £ = 0, it 
becomes Xmin = 3 - 6 f° r £ = 4% and Xmin = 9 for ^ =8%. 

In conclusion we have shown that the MSW solution to the SNP exists for any realistic levels of 
matter density noise (£ < 4%). Moreover the MSW solution is essentially stable in mass (4-10 _6 eV 2 < 
5m 2 < 10" 5 eV 2 at 90% CL), whereas the mixing appears more sensitive to the level of fluctuations. 



6. We can reverse our point of view, wondering whether the solar neutrino experiments can be a 
tool to get information on the the level of matter noise in the sun. In particular, the future Borexino 



experiment [15], aiming to detect the 7 Be neutrino flux, could be sensitive to the presence of solar 
matter fluctuations. In the relevant MSW parameter region for the noiseless case, the Borexino signal 
cannot be definitely predicted (see Fig. 3a). Within the present allowed C.L. regions (dotted line) 
the expected rate, Z Be = R V bT I 'R-Be^ ( son d lines), is in the range 0.2 -j- 0.7. 

On the other hand, when the matter density noise is switched on, e.g. ^ = 4% (see Fig. 3b), 
the minimal allowed value for Zs e becomes higher, Zs e > 0.4. Hence, if the MSW mechanism is 
responsible for the solar neutrino deficit and Borexino experiment detects a low signal, say Zs e ^ 0.3 
(with good accuracy) this will imply that a 4% level of matter fluctuations in the central region of 
the sun is unlikely. The same argument can be applied to v e — > u s resonant conversion, whenever 
future large detectors such as Super-Kamiokande and/or the Sudbury Neutrino Observatory (SNO) 
establish through, e.g. the measurement of the charged to neutral current ratio, that the deficit 
of solar neutrinos is due to this kind of transition. The expected signal in Borexino is very small 
Zse ~ 0.02 for £ = (see Fig. 3c). On the other hand with £ = 4%, the minimum expected Borexino 
signal is 10 times higher than in the noiseless case, so that if Borexino detects a rate Zse ^0.1 (see 
Fig. 3d) this would again exclude noise levels above 4%. 

Let us notice that Super-Kamiokande and SNO experiments, being sensitive only to the higher 
energy Boron neutrinos, probably do not offer similar possibility to probe such matter fluctuations 
in the sun. 

The previous discussion, which certainly deserves a more accurate analysis involving also the 
theoretical uncertainties in the 7 Be neutrino flux, shows the close link between neutrino physics and 
solar physics. 

This work has been supported by the grant N. ERBCHBI CT-941592 of the Human Capital and 
Mobility Program. 
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Fig. 1: The averaged solar neutrino survival probability P versus E/5m 2 for small mixing angle, 
sin 2 29 = 0.01, (Fig. la) and for large mixing angle, sin 2 29 = 0.7, (Fig. lb). The different curves 
refer to different values of matter noise level £ as indicated. 
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Fig. 2: The 90% C.L. allowed regions for the v e — ► v^ T (Fig. 2a) and for the v e — > v s (Fig. 2b) 
conversion. The different curves refer to different values of matter noise level £ as indicated. 
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Fig. 3: The iso Zs e = R P Be /^Be 95 contours (figures at curve) in the v — e scattering Borexino 
detector (solid lines). The threshold energy for the recoil electron detection is 0.25 MeV. The 90% 
C.L. regions (dotted line) and the corresponding best fit point are also drawn. Fig. 3a and Fig. 3b 
refer to the case of v e — > v^ T conversion and for £ = and £ = 4%, respectively. Fig. 3c and Fig. 3d 
refer to the case of v e — ► v s conversion and for £ = and £ = 4%, respectively. 



